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We discuss how measurements of fluctuations in the absorption of cosmic microwave background 
(CMB) photons by neutral gas during the cosmic dark ages, at redshifts z « 7-200, could re- 
veal the primordial deuterium abundance of the Universe. The strength of the cross-correlation of 
brightness-temperature fluctuations due to resonant absorption of CMB photons in the 21-cm line of 
neutral hydrogen with those due to resonant absorption of CMB photons in the 92-cm line of neutral 
deuterium is proportional to the fossil deuterium to hydrogen ratio [D/H] fixed during big bang nu- 
cleosynthesis (BBN). Although technically challenging, this measurement could provide the cleanest 
possible determination of [D/H], free from contamination by structure formation processes at lower 
redshifts, and has the potential to improve BBN constraints to the baryon density of the Universe 
Sib ft 2 - We also present our results for the thermal spin-change cross-section for deuterium-hydrogen 
scattering, which may be useful in a more general context than we describe here. 
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Introduction — After the cosmic microwave background 
(CMB) radiation decoupled from the baryons at a red- 
shift z w 1100, most CMB photons propagated unfet- 
tered through the neutral primordial medium. This has 
allowed exquisite measurements of the temperature fluc- 
tuations in the primordial plasma at the surface of last 
scattering, and the statistical properties of these fluctua- 
tions have recently been used, in conjunction with other 
observations, to determine the cosmology of our Uni- 
verse 0. After the photons kinetically decoupled from 
the gas at z ~ 200, the latter cooled adiabatically with 
T g oc (1 + z) 2 , faster than the T 7 oc (1 + z) cooling of the 
CMB. This epoch, with most of the baryons in the form 
of relatively cold neutral atoms and before the first stars 
formed, is known as the cosmic dark ages. 

The reason most CMB photons propagate unimpeded 
through the neutral primordial gas is elementary quan- 
tum mechanics — atoms absorb non-ionizing radiation 
only at the discrete wavelengths determined by the dif- 
ferences of their atomic energy levels. One interesting 
example is the well-known 21-cm spin- flip transition [2|, 
due to the hyperfine splitting of the ground state of the 
hydrogen (H) atom. At any given z, CMB photons with 
wavelength A21 = 21.1 cm can resonantly excite this 
transition. By measuring brightness-temperature fluctu- 
ations due to density fluctuations in the neutral gas [3j, 
radio telescopes observing at A = (1 + z)Xqi can probe 
the matter power spectrum at z » 30-200 [4J. 

In this Letter we discuss another application of these 
measurements. Less well-known than the 21-cm transi- 
tion of neutral H is the spin-flip transition of neutral deu- 
terium (D) at A92 = 91.6 cm 0,1^. We show below that 
cross-correlating brightness-temperature fluctuations at 
a wavelength A H = (1 + z)\ 2 i with those at a wavelength 
A D = (l+z)Ag2 allows a measurement of the primordial D 
abundance. In principle, this technique could constrain 
the primordial value of [D/H] = n D /n H to better than 
1%. While there is no physical obstacle to such a mea- 



surement, it would certainly be technically challenging, 
and require a heroic experimental effort; simply detect- 
ing neutral D during the cosmic dark ages would be a 
significantly easier goal. 

Deuterium has long been recognized as our best 'bary- 
ometer' because its primeval relic abundance is so sensi- 
tive to the baryon-to-photon ratio jn = n^/s. Moreover, 
big bang nucleosynthesis (BBN) |jj is the only known 
natural production mechanism, although mechanisms in- 
side galaxies can destroy it || . The measurement we de- 
scribe below could thus determine the true BBN abun- 
dance of D and, in principle, might improve BBN con- 
straints to the baryon density of the Universe fl^h 2 . 
Hyperfine Structure of H and D Atoms — The \i ■ B in- 
teraction between the magnetic moments of the electron 
and the nucleus splits the ground state of single-electron 
atoms into eigenstates of the total spin operator F = S+I 
with eigenvalues F + = 1+1/2 and F_ = 1 — 1/2 and 
AE = (16/3)F + M. B ( 5N M. N /ag) (e.g., 0). Here, S is elec- 
tron spin, I is nuclear spin, ao is the Bohr radius, u. B is 
the Bohr magneton, [i N is the nuclear magneton, and <? N 
is the nuclear g factor (g p = 5.56 for H; g D = 0.857 for 
D). The proton, with / = 1/2, splits the H ground state 
into a triplet with F + = 1 and a singlet with F_ = 0. 
The deuteron, with 7 = 1, splits the D ground state into 
a quartet with F + — 3/2 and a doublet with F_ = 1/2. 

The Spin- Temperature — The population of atoms in the 
excited spin state relative to the ground state n + /n_ = 
Q7+/<7_)exp{— T^/Tg) can be characterized by a spin tem- 
perature T s . Here g + = 2F + + 1 and g_ = 2F_ + 1 are 
the spin degeneracy factors and T* = AE/k B . For H and 
D, we have respectively T? = 0.0682 K, T° = 0.0157 K, 
and (g«/g«) = 3, (g?/ 5 D ) = 2. 

Competition between three factors determines T s : ab- 
sorption of 21-cm CMB photons, absorption and re- 
emission of Lyman-a photons (the Wouthuysen-Field or 
WF effect 0, EI), and atomic spin-change collisions 
(collisions with free electrons are unimportant in these 
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environments |Icj). The first drives T s toward T 7 , while 
the latter two drive it toward the gas temperature T g . 
In equilibrium the spin temperature of the a species X 
(either D or H) is 



(l + X x )T g T 7 
(T g + X X T 7 ) ' 



(1) 



where x x = Xc + Xa is the sum of the equilibrium 
threshold parameters for spin-change collisions and for 
radiative coupling through the WF effect. Explicitly, 
X? - (Cf_T?)/(A*_T 7 ) and xl = {P+- T x ) / (A*_ T 7 ) , 
where C*_ is the collisional de-excitation rate, A*_ is an 
Einstein coefficient, and Pf_ oc P a , where P a is the to- 
tal Lyman-a scattering rate. At z ^> 10, before the first 
galaxies formed, P a is tiny and the WF effect can be ne- 
glected. However, it might have interesting consequences 
near z ~ 10. 

H-H and D-H Collision Rates — While the cross section 
for H- H sp in-change collisions cr™ is well known |12l Il3l 
0- E 113 ■ we were unable to locate the D-H spin-change 
cross section for the temperature range of interest and 
computed cr°" using standard methods 1 . In the elastic 
approximation, 
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FIG. 1: The thermal spin-change cross sections that keep 
T s collisionally coupled to T 9 for D (solid line), H (dashed 
line), and muonium (dotted line). Although the potentials 
are identical, the peaks differ because of the reduced masses. 
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where k = fj, m v/h, /i DH = m D m H /(m D + m H ) is the re- 
duced mass of the D-H system, and v is the relative ve- 
locity. 2 The partial wave phase shifts in the triplet and 
singlet electronic potentials V t (R) and V S (R) in which 
the D and H atoms scatter are i rjf H and s r/f H respec- 
tively. We used the variational potentials of Refs. [2lll22l | 
for R < 12 a , the smooth fit of Ref. 23] for V t (R) 
between 12 ao < R < 15 cto, and the asymptotic form 
of Ref. [24| for larger R. We found the phase shifts 
by solving the radial Schrodinger equation at each en- 
ergy E = h 2 k 2 /2/i DH and angular momentum y/l(l + 1)H, 
truncating at l ma x ~ 100 to resolve the resonant struc- 



ture at sufficiently high E. We verified that our results 

unting for our n 
recent potentials), with those of Refs. 0, 
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agree with Ref. [19| at 1 K (after accounting for our more 
;hose of Refs. [3, EH EE E3 3 for 
for the muonium-H spin-change 



ct™, and with Ref. 1 
cross section cr^ H . 



The collisional de-excitation rate is C* = % H (T* 



where v XH = ^J^k^Tgj '/j XH 7r is the thermal velocity, a 
is the thermal spin-change cross-section (averaged over 
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the Maxwell-Boltzmann distribution of relative veloci- 
ties), and n H is the number density of H atoms. In 
Fig. fflwe plot a-™, 5-° H , and e^ H . While a™ falls off 
for T g < 100 K, continues to rise to a peak near 
T g ~ 1 K. This occurs because of low-energy s-wave and 
p-wave contributions to D-H scattering (a much larger 
scattering length and a scattering resonance). These do 
not appear for H-H because of the differing reduced mass 
(Mhh ~ while /j. D h_~ 2to h /3). The discussion of D- 

H spin-change in Ref. [lfj did not account for this and 
incorrectly concluded that cr°" ~ a™. 
Spin- Temperature Evolution In Fig. |2]we plot T 7 , T g 
(found using reef ast HH), T S H , and T S D as a function of 
z. After the gas cools below T 7 , collisions keep T S H and 
T S D coupled to T g . Near z ~ 30 collisions become ineffi- 
cient for H and T S H returns to T 7 . T S D remains coupled 
to T g down to significantly lower redshift both because 
the lifetime of the excited state of D is relatively long 
(Al_/A°_ = 61.35) and because a™ > ct™ at low tem- 
peratures. 

Brightness Temperature Fluctuations — When the spin 
temperature of a given species is less than T 7 it will 
absorb CMB photons. The brightness temperature is 
T b x = ar x (T s x - T 7 ), where 



1 It has b een computed at higher temperatur es b y Ref. llSl , at 1 K 
by Ref. llSl . and measured at 1 K by Ref. |2(1 . 

2 Compare with cr™ = (7r/4fc 2 ) yjg (2« + l)sin 2 (^f H - a r;f h) for 
H-H collisions when quantum sym metry can be neglected llll . 

3 We are in harmony with Ref. 1 1(1 which found that — ► 0) = 
4/3 larger than previously quoted Il.'lll4l at high T g . 
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is the optical depth of the spin-flip transition in ques- 
tion, a = 1/(1 + z) and TC(z) is the Hubble parame- 
ter. We are interested in correlations between brightness 
temperature fluctuations S^, (n, a) = ST£(h,a)/T£(a) = 
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FIG. 2: The H and D spin temperatures as a function of z. 
Here we assume P a — for all z. 



Pj, (a)S(n,a) observed in a direction n at wavelengths 
differing by a factor A92/A21. Here, 
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<5(n, a) = <5n H (n, a)/n H (a) = <5n D (n, a)/n D (a) is the den- 
sity contrast, and % x = x x (l + X x )- At high z, when 
T g ss T 7 , r — > due to residual Thomson scattering with 
free electrons [2(| (fluctuations are isothermal), but as 
the gas begins to cool adiabatically T — > 2/3 |27|. In 
Eq. Q we have neglected the contributions to 5r b from 
fluctuations in the neutral fraction (likely to be small at 
high z) and, for simplicity, fluctuations in the gradient of 
the radial velocity Sg rVr |2jg . The latter will enhance our 
signal by a factor of ~ 1-2. In Fig. El we plot T b , a/?f b T£, 
f b D and a/?£ b f b d , where T b D = eT b D and e = [D/H]. We see 
that and T b peak at much lower z than their H 
counterparts because, as discussed above, T S D is coupled 
to T g to lower z. 

D-H Cross Correlations — We now estimate the cross- 
correlation of brightness temperature fluctuations across 
frequencies related by A21/A92. We write the brightness 
temperature fluctuation due to H or D as H(h, a) — 
(3^ b (a)T^(a)S(h,a) and eD(h,a) = e(3» b (a)fg(a)S(h,a) 
respectively. A radio telescope observing at a frequency 
v will measure the quantity C[n;^] = H(h, v/v2i) + 
e£>(n, v/vq2) + N[ii; v], where N[h; v] is the instrumen- 
tal noise. We form the product C[n; v a ]0[h] r/g], where 
Vp = {v<jijv-2,\)v a . Assuming that <$(n, a) is a zero-mean 
Gaussian random field and uncorrelated Gaussian noise, 
its expectation value is (0[h; v a ]0{h; vp\) = e(H a Df}) 
to leading order in e. Here we have introduced the 
shorthand H a = H(n, v a /v2i), Dp = -D(n, vp/v^), and 
N a = N[h-v a \. 
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FIG. 3: The brightness temperatures Tjf and T° and the 
brightness temperature fluctuations /3^ h T^ and f3° h T° (scaled 
to the growth rate of density perturbations 5 oc a). 



We now understand the crucial point of this Letter. 
The 21-cm and 92- cm fluctuations at these frequency sep- 
arations must be correlated because they trace the same 
underlying patches of the Universe. 

Note that we have neglected the relatively small in- 
trinsic correlations of the H brightness fluctuations from 
large-scale modes of the density field which contribute at 
the level of ~0.1% or less of the D signal. If necessary, 
these intrinsic H correlations could be removed through 
independent measurements of the matter power spectrum 
or by correlating the 21-cm signal with frequencies near 
but not equal to the corresponding D patches. 
Signal Estimate — The signal-to-noise contributed by a 
pair of frequency bands centered around {v a ,vp) for an 
experiment with a maximum baseline of L and frequency 
resolution Ai> is 
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(5) 



where Op = A/3 / L is the angular resolution in the D 
band. Here, (H a D p ) = <J 2 S ((3^)(^T°) and (H 2 a ) = 
^(0fjbf, where 



dkkP(k) / dk x jfak z p) 1 v 
Jo 



kip) 



(fc 2 - kl)p 2 



(6) 



is the variance in the density field at scale-factor a 
smoothed over the coin-shaped regions of the Universe 
of comoving radius p and thickness 2£ p that arc sampled 
by an experiment with high spectral (radial) resolution 
but lower angular (transverse) resolution. We adopt a 

T s 2 ys /(/ c 2 ov AW int ) (e.g., 0), 



noise variance of (A~ 2 ) 
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where T sys = 6500[^ Q /(30MHz)]" 2 K is the noise tem- 
perature, / cov is the covering fraction of the array, and 
tint is the integration time. Our choice for T sys is only an 
estimate, and the noise (ultimately due to Galactic syn- 
chrotron radiation) varies strongly across the sky. The 
total signal-to- noise ratio (S/M) to t is the sum in quadra- 
ture over all pairs of frequency bands (f a ,u^). 

If collisions dominate the coupling between T S H and T g 
down to z ~ 7 then a value [D/H] ~ 3 x 10~ 5 could be 
detected at 1- to 2-a by an experiment with L ~ 7.5 km 
and Av ~ 100 kHz in ~6 years. If, however, the first 
generation of stars created a flux of Lyman-a photons 
which coupled T S H to T g until z ~ 7 through the WF effect 
without significantly heating the g similar detection 
might be made by a smaller experiment with L ~ 2.5 km 
and Av ~ 1 kHz. Although it dramatically enhances the 
21-cm signal, the WF effect does not improve the [D/H] 
measurement by the same margin because S/Af becomes 
independent of H a once (H^) ^> (N%) (it only serves 
to make the H fluctuations a better matched template). 
Finally, we note that an experiment capable of mapping 
21-cm brightness-temperature fluctuations out to l rnax ~ 
10 5 (where it may be a powerful probe of the small-scale 
matter power spectrum Q) could measure [D/H] to a 
precision as good as ~ 1% — or even ~ 0.1% if the WF 
effect coupling is efficient. 

For these estimates we have assumed a ACDM cos- 
mology with n s = 1, and that a significant fraction of the 
Universe remains neutral until z ~ 7. The largest contri- 
bution to the signal originates from z < 10 where the D 
signal peaks and the variance in the density fluctuations 
is largest. Varying these assumptions, or including addi- 
tional sources of temperature fluctuations, could change 
(S/Af) tot by factors of order unity. 

Discussion — Despite the obvious technical challenges in 
observing this signal, we emphasize that it has the virtue 
of providing the cleanest possible measurement of the 
primordial [D/H], free from contamination by structure 
formation processes at lower z. Via the window of BBN, 
this would allow radio telescopes to peer into the first few 
minutes of the Universe. We believe future searches for 
cosmic 21-cm fluctuations should bear this possibility in 
mind. 

We also note that 3 He + has a hyperfine transition 
(with A = 3.46 cm) that can be used in a similar fashion; 
it has the advantage of much lower foreground contami- 
nation at higher frequencies. This line will appear during 
reionization and should exhibit a strong anti-correlation 
with the corresponding 21-cm signal. If the astrophysics 
of reionization can be understood well enough, the cross- 
correlation of this line with the 21-cm line could supple- 
ment the D-H experiment in order to probe BBN in even 
more detail. 
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